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We look at the possibihty of non-perturbative particle production after inflation from SUSY flat 
directions produced by rotating eigenstates thereby avoiding the standard adiabaticity conditions. 
This might lead to preheating and prevent the delay of thermalisation of the universe. We investigate 
the flat directions LLE'^ and U'^D'^D'^ and find no particle production. These 2 directions are very 
important, since they have been named as possible candidates for being the inflaton. We investigate 
QLQLQLE'^ and find particle production and therefore the possibility of preheating. We investigate 
the LLE'^ and U'^D'^D'^ directions appearing simultaneously, and find no production. Finally, we 
investigate LLE'^ and QLD'^ simultaneously - with one L-field in common. Here we do find particle 
production and therefore the possibility of preheating. This means that if SUSY fiat directions are 
to delay thermalisation and thus explain the (lack of) gravitino production, it is necessary to explain 
why complicated directions as QLQLQLE'^ are not exited, and why combinations like LLE'^ and 
QLD'^ axe not both exited. 
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I. INTRODUCTION 



The scalar potential of the Minimal Supersymmetric Standard Model (MSSM) possesses a large number of F- and 
D-flat directions P, 0. These flat directions might generate the baryon asymmetry of the Universe through the 
out-of-equilibrium CP violating decay of coherent field oscillations along the flat directions themselves ■ 

Recently the cosmological importance of flat direction vacuum expectation values (VEV)s [Often we will use VEV 
meaning nonzero vacuum expectation value - this should not cause confusion] and the decay thereof has been inves- 
tigated. In ^ it was asserted that large fiat direction VEV's can persist long enough to delay thermalization after 
inflation and therefore lead to low reheat temperatures. This is of great importance. A lower reheating temperature 
would potentially solve the (lack of) gravitino problemQ. It has also been claimed Q that large flat direction VEV's 
can prevent non-perturbative parametric resonant decay (preheating) of the inflaton since the inflaton decay products 
become sufficiently massive preventing preheating from ever becoming efficient. These arguments hold so long as the 
fiat direction VEV's do not rapidly decay - they must persist long enough so that they can delay thermalization and 
block infiaton preheating. In [9| it was claimed that non-perturbative decay can lead to a rapid depletion of the fiat 
direction condensate and thus precludes the delay of thermalization after infiation. It was also concluded that in 
order for the flat direction to decay non-perturbatively the system requires more than one flat direction 1,13. In 
[1^ it was stated that even in the presence of multiple flat directions, some degree of flne-tuning was necessary to 
achieve flat direction decay. We note that very recently [ll| it has been claimed that even if non-perturbative particle 
production happen, the main decay mode will still be perturbative. 

The presence of Nambu-Goldstone bosons (Goldstones) is a very important point in this discussion. The fiat 
directions are charged under the gauge group of the MSSM. Therefore the fiat direction VEV will break some or all 
of the gauge symmetries of the theory and therefore the presence of the associated Goldstones must be expected. 
[qI considers a gauged C/(l) model and constructs the mixing matrix for the excitations around the fiat direction 
VEV. In it was claimed that in the single flat direction case, non-perturbative decay proceeds solely via a massless 
Goldstone mode as only the Goldstone mode mixes with the Higgs and all other massless moduli remain decoupled. 
Since the Goldstone represents an unphysical gauge degree of freedom, it was concluded [1, that no preheating 
occurs in the single fiat direction case - a Goldstone can be gauged away. In order to determine if flat direction VEV's 
decay non-perturbatively into scalar degrees of freedom, one must remove the Goldstones ie. use the unitary gauge. 

In an earlier paper we and our colleagues considered toy models to demonstrate that, in the unitary gauge, 
the mixing matrix of the excitations around a fiat direction VEV permits preheating. Moreover, we found that fiat 
direction decay depends on the number of dynamical, physical phases appearing in the flat direction VEV. Speciflcally, 
a physical phase difference between two of the individual fleld VEV's making up the flat direction is needed. 

In the present paper, we look at (some of) the actual SUSY flat directions. In section |ll] we look at the LLE^ 
direction, in section IIIII we review the particle production mechanism from rotating eigenvectors, in section IIVI we 
conclude on the LLE'^ case, and in section|V]we do the U'^D'^D'^. The mentioned directions are especially interesting, 
since they are mentioned as especially well suited infiaton candidates in Q. We investigate QLQLQLE^ in section 
IVII We conclude on one flat direction in section IVIII Then we proceed to 2 directions, first the non-overlapping 
LLE'^ + U'^D'^D'^ of the 2 inflaton candidates of [6] in section IVIIII and then the overlapping directions QLD'^ + LLE'^ 
in section |Xl Finally we look at the simpler approach of just counting the flclds without any calculations in section 
IXII and conclude in section IXIII 



II. LLE" 

One flat direction often mentioned in the literature is LLE'^. Flatness demands the fields with VEV's to come from 
different generations, and the 2 L fields with VEV to have opposite SU(2)-charge. Also, the 3 VEV's must have the 
same absolute value. This leaves essentially only 1 choice (when masses are ignored). 

We give these VEV's: 

<fi>= ipe"''' (1) 

Also, it is clear that 2 other fields will play a role. 



< e >= 

; >= 0- (2) 



The Lagrangian reads 
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C = Y l\DM' -V- - jK^ (3) 



1=1 
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where for field (bi 



bj denotes the covariant derivative. P° is the a*'* Pauli- 



matrix. The potential we consider arises from the supersymmetric D-terms and has the form 
where 

i 

Da - |0^F'^</> (6) 

where qi is the hypercharge, and 51,52 are the hypercharge and SU2 gauge couplings. 

The essential part is removing the Goldstones correctly. To do that we start by looking at the fields with the VEV's 
only (no excitations). We've written those earlier, and we get mixed kinetic terms 

C D -ip^Ao{<7l + (72 - 2(73) - y^^A3{<7l - (72) (7) 

which has the form of a coupling between the gauge field and the background condensate. Terms of this type will 
feed into the equations of motion for the gauge field which, in turn, will have an effect on the equations of motion for 
the scalar excitations. The remaining terms in £ are 

^ifi^ [6Al + 2A\ + + 2Al + gx^ + 62" + (73^] (8) 
- all desired terms. By making a ?7(1) gauge transformation on the VEV, 

^ {^\) = e'"-^ (9) 

with 

A=^^^^f^, (10) 
and by making a SU{2) gauge transformation on the VEV, 

^ ($:) = e'^'^ (11) 

with 

7=^, (12) 

we can gauge the unwanted terms away and avoid a complicated analysis of the kinetic terms. The resulting form of 
the VEV reads, 

(m) = ipe^" (13) 

where a — {ai + (72 + (73)/3 represents the remaining independent physical phase. Following [l3l |. we can write the 
fields in the unitary gauge as (including the other relevant fields), 

j^e = (^ + 6)e'^"+^' 

- {^^ + ^£,&)e'" (14) 

= (^ + e3)e^^"+*^ 

^ ((^ + ^4)6'*"^^^ 

where a represents time dependent phase of the VEV (we have just showed the phase differences are gauged away), 
5i parameterises its excitation, ^2,3,4 parameterise the excitations around the VEV, and ^5,5,7,8 parameterise the 2 
no- VEV fields (the phase on these fields is not necessary, but allowed, and will be convenient). 
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Again we will look at the kinetic term. First, the (/j^-term 

[6^2 + 2Al + 2Al + 2Al + Sct^] (15) 

- not surprisingly. This contains no goldstones, so we proceed to next order. 
The terms indicating Goldstones should include ^i. These terms are 

£3-^(^1(^6 + 4) + ^2(^7-^5)). (16) 

The remaining terms are on the forms 

£ D (c^jkirAjAk + D.j^.Ajd + E.^.a"^ - Vs^i.tr) (17) 

Here it is clear, that the field excitation terms (excluding derivative terms) are suppressed compared to iy9^-terms in 
the potential. The only excitation is the last term. However, it is just a "mixing" between an excitation and its own 
VEV. 

The Goldstones are removed by demanding 



^8 and ^7 = ^5. Doing this, and renormalising, we take 

(^ + ^2)6^^^+^^ 
(Cs + i^e) 

(gs - ij^) 



(18) 



This does indeed kill the mixed derivative terms. The remaining terms stay as they are. But they are all VEV- 
suppressed, so it is justified to move to the coordinate derivative, rather than the covariant derivative. 
The remaining kinetic term (to zero'th order in ip) are 



i=i 



4 

■E 

1=2 



(19) 



The first term are the kinetic terms that show we have correctly normalised kinetic fields - including that there are no 
cross terms. The second term is completely negligible compared to the (p^-terms of V. The third term, though bigger 
than the prior one, is still suppressed. The fourth term is a rotation between the excitation states. These are very 
important and will give the U- matrix below. The fifth term is just a VEV-term, and the final term is the "mixing" 
between an excitation and its own VEV. So everything is fine. 

On substituting the fields of eqUHinto the Lagrangian given in eq|3]and defining the vector S = (^1,^2,^3: Csi Ce)'^; 
we find the quadratic terms 



(20) 



where the ellipses denote higher order terms and interactions. The matrix U given in the second part of the third 
term in eg 1 191 reads 



(0 


2(T 

V3 


2(T 

V3 


2(T 

V3 





o\ 
























































































0.) 



(21) 



However, we want an antisymmetric matrix for the procedure below. Using partial integration - and ignoring surface 
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terms - we find 



U : 



/ --2- --2- --2- 
/ " .AJ .AJ 








& 

V3 


& 

V3 
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a 

V3 






































V3 




























































(22) 



while the mass matrix for the physical excitations appears as 
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9l 


+ 9l 


9l 


-9l 


-25? 











9l 


-9l 


9l 


+ 9l 


-25? 















^^91 






45? 


























25i 







Vo 


















25i 


/ 



bmIb'^, 



(23) 



with eigenvalues — &gi^p , Af| = M| = Af| = 2glip^,M^ = Mg = (the entries of the diagonal matrix Aid)- B is 
an orthogonal matrix which diagonalises Ai^ and M^_4 corresponds to the mass of the physical eigenstates associated 
with the spontaneous breaking of the symmetries. M| = Af| = correspond to the massless excitations around the 
fiat direction VEV. 

The last term in eg 1201 appears as a consequence of the time-dependence of the background ~ it represents a mixing 
between the fields ^i, 2, 3,4, 5, 6) and their time-derivatives. The effect of these terms on the system becomes clear if we 
make field redefinitions that remove the mixed derivative terms. The resulting transformation leaves the system in 
an inertial frame in field space and leads to a time-dependent mass matrix. Defining S' = AS (A is orthogonal), we 
find the condition that A must satisfy in order for all the mixed derivative terms to cancel 



A^A = U. (24) 

The Lagrangian for the S' system now reads 

C^\\d,^r-\^'^M''E' (25) 

where M''^ = AM'^A^ = ABM'^B'^ = CM'^C'^, and C = AB. The matrix C is an orthogonal time-dependent 
matrix, with columns corresponding to the eigenvectors of A^'^. We now have a system of scalar fields with canonically 
normalized kinetic terms and time dependent eigenvectors. 

The central point of this discussion centers precisely on the appearance of the time dependent eigenvectors for the 
six scalar fields. This satisfies a necessary but not sufficient condition for preheating. In the next section, we briefly 
run through the details of the non-pert urbative production of the light scalar fields following the analysis of il4]. This 
is a more brief summary than in [12'] - which is otherwise followed here. 



III. NON-PERTURBATIVE PRODUCTION OF PARTICLES 



Including gravity, the dynamics of the re-scaled conformally coupled scalar fields, Xi = a^J, where a denotes the 
scale factor and the i-th component of the vector S', are governed by the following equations of motion (sum over 
repeated indices is implied), 

x^ + nl{t)xJ^o (26) 

where dots represent derivatives with respect to conformal time t, and 

nl=a^M'l + k^6,„ (27) 

where k labels the comoving momentum. Using an orthogonal time-dependent matrix C(t), we can diagonalise fiy 
via C'^(t)ri^(t)C(i) = (j-i^it), giving the diagonal entries ujj{t). Terms of the form ^ fi^X arising from the kinetic 
terms do not affect the evolution of the nonzero k quantum modes [l^ . 

As the vacuum changes, a new set of creation/annihilation operators are required. We use Bogolyubov transforma- 
tion with Bogolyubov coefficients a and /? (which denote matrices in the multi- field case). 
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Initially a = I and /3 — while the coupled differential equations (matrix multiplication implied): 



a = ^iLoa H (3 — la — J fi 

P = -^a + iuj/3 - Ja- 1/3, 

ZUJ 

govern the system's time evolution with the matrices / and J given by 

2 \ Juj Ju) 



(28) 



(29) 



(30) 



Similarly to the single-field case it can be shown that at any generic time the occupation number of the ith 
bosonic eigenstate reads (no summation implied) 



(31) 



As pointed out in [13], there exists two sources of non-adiabaticity in the multi-field scenario. The first source 
arises from the individual frequency time dependence and appears as the only source of non-adiabaticity in the single 
field case. The second source appears from the time dependence of the frequency matrix J7y giving rise to terms in 
eg 1281 proportional to / and J. This second source provides the most important contribution in our analysis and gives 
rise to non-perturbative particle production. 

Since initially a = I and /? = 0, eg 1281 shows that a non- vanishing matrix J is a necessary condition to obtain /3 7^ 
and hence ni{t) 7^ 0. In general, we have 



C'^C = B^A^Ab = -B^UB 



(32) 



where A, B and U were defined in the previous section. The last equation only holds if B is constant in time. This 
is obviously the case in the LLE'^-case, since, B diagonalises a constant matrix. 

IV. LLE" CONCLUSION 

For the LLE'^ example outlined above, J is a 6 x 6 zero matrix. Therefore there is no particle production and no 
preheating. 

V. VD^D" 

One would expect the U'^D'^D'^ case to be much the same - as indeed we shall see it is. We give VEV's to these 
fields (again from different generations to avoid F-terms) 



(33) 



The Lagrangian reads 



(34) 



where for field 4^ 
A^^ Gell-Mann-matrix. The potential now looks like 



4 ^ 4 - 

i 

denotes the covariant derivative, where GM^ is the 



(35) 
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where 



Da = ^cji^GM^^ 



(36) 



where 53 is SU(3) gauge coupUngs. Removing mixed kinetic terms as before (in 2 tempi), we use 



u^i = (^ + ^2)e'^''+^\ 

(?5 + iS.fi) Aa 



,c2 



V2 



^ (ii+zes)^, 

\/2 



„c2 



V2 



gc3 ^ (6+iil£)gio 
.cT _ (^7 - ia- 



V2 



(37) 



The U-matrix is 



U = 



( 





dr 

V3 


(j 

V3 


a 

V3 

















0\ 
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V3 
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v{3 































(T 





























































































































































































\ 
































(38) 



while the mass matrix for the physical excitations appears as 
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3 53 
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253V 



BMlB^, (39) 



with eigenvalues = (§52 + 2gl) ip"^, M| = M| = M| = M| = M| = M| = M| = 2glip^,M^ = Mf^ = 0. Also 
here, we end up with J = and no particle production and therefore preheating. 



VI. QLQLQLE" 



QLQLQLE" case has so many fields with nonzero VEV, that all the phase differences cannot be gauged away. The 
starting point could be (notice that here, there are 2 essentially different possibilities - the squarks having identical 
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SU (2)-charge - or not) 



The Lagrangian reads 



< u^^ >= 


(fie 


< 0^=2 > = 


(fie 


< 4^=^ >= 


(fie 


< e >= 


(fie 


< /U > = 


(fie 


< T > = 


(fie 


<e^ >= 


(fie 



i<72 
i<73 



(40) 



where for field <pi = 
potential now looks like 



i=l i 



(41) 



denotes the covariant derivative. The 



(42) 
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; terms as before (in 2 tempi), we use 





(^ + ^4)6*^"^ + ^' 








d'^ = 


(^ + ^5)6^^"^^^' 








t"^ = 










e = 


(V + C7)e^^'^^ + ^^ 


-cr2 + 






M = 


(^ + ^8)^^'^"^ + ^" 


-<72- 


^+^^ + 


_S3_) 


T = 


(<p + 6)e*('^^ + ^" 


-2(73 


V649 ' 




— 




) 






= 


x/2 








U^^ = 


V2 









^cl ^ 65 + i^l6 _ 61 + ^62 _ C27 + »68 _^ 69 - ^60 _^ 61 " ^$32 \ giCTi 

V 72 76 2V3 2^5 V30 J 



^c3 ^ 69 + »C20 ^»ai 
72 



67 + ^68 C29 - ^60 61 - ^62 \ 




V 2^1 J 




I I 69 + ^$30 _ 61 + ^62 \ gi(o-i-2o-3) 

V a/2 76 2V3 2^5 V30 / 

The non-zero elements of the U-matrix are 



Uai = U51 — Ue,i — UioA = —7= 



U21 



22 



U23. 



24 



C^l - C^2 +<^3 

V7 
a/2 

- ^2 + 

a/7 
crj - 20-3 



t^7,2 
«77,3 

U30 



(^20.15 — 



V3 
0^1 - 0^3 

ai - f73 



20,15 = i^iejig = 

t^22,15 = t^l6,21 = 
C^24,15 = C^23,16 = 

V 10 

1^26,15 — 1^25,16 
t^l9,20 



'lb 

2a[ + a'3 



U22A9 = f^20,21 



3 

CTi — (To 



3^2 

t^24,19 = [^23,20 = /-- 

V o(J 

1^26,19 — 1^25,20 



3^5 
5(7^ + (73 



C^24,21 = f^23,22 = 



2V15 

1 — 3 

t^26,21 = C^25,22 ^ 



C^26,23 — [^24,25 
t^25,26 



3V10 

7(7^ + 8(7^ + 7(7^ 

10 

3<7i + 2(7^ + 3(7^j 

5\/6 

4cri - 4(72 + 4(73 



and their antisymmetric counterparts. 
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The mass matrix for the physical excitations appears as (in units of 

9 3 



M4,4 


= M5,5 = 


M6,6 




= M4,6 = 






= M4,8 = 


M4,9 


M4,io 


= M5,io = 


3 


Mr,7 


= M7,8 = 


M7,9 



M7,10 


= =M8,10 = 


= M9,io = 


Mio.io 


= 4<7? 






= Mi2,12 = 






= Mi6,l6 = 


252 




= Afl6,20 = 


25i 
^/3 



ffl , 2 25I 

M5,7 = M5,8 = M5,9 = M6,7 = M6,8 = ^5,9 = ^ - §1 



Ms,s = M8,9 = M9,9 = gl + gi 

„2 



-^15,21 = Afl6,22 = 



-^15, 25 — —-^16^26 



^/3 
3\^ 

2\/3.9i 



Mi5,23 = -Mi6,24 = (45) 



-^19,19 = Af20,20 = -P- 



3 

Afi9 21 = M20,22 = 



3 



-^19,23 = —-^20,24 = 

-^19,25 = — Af20,26 = 

-^21,21 = Af22,22 = ^ 

-^21,23 = — M22,24 = ^^^^ 

M - M - 

-'«21,25 — —■''^22,26 — ^ 



95,2 

-^23,23 = -^24,24 ^ 



-^23,25 = M24,26 



2 

5 

3\/65i 



-^25,25 — -^26,26 



5 

and their symmetric counterparts. The eigenvalues arc Mf = ^'^^i+Sgs+y/i'^ig^-Sig^g^+sig^ ^2^ ^2 _ 

llgl+9gl-^121gt-5igfgi+Slgt ^2^M2 ^ ^2 ^ GffaV^^ M| = M| = = Mi = = M^g = Mf^ = M^^ = 

2glv\Mis = .... = Mie = 0. 

The J-matrix is (really: The J-matrix can be the splitting of eigenspaces of higher dimensions into subspaces is 
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arbitrary) 



16 = ~Ja 15 = — ^ r (ScTo + GCTq) 



^4,17 = ■/a, 
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-Vk+\ k 



''''' - - ^4V2(P + 652V) 



9fc^+3(llg;+9g|+^121gf-54g;g|+81^)y'= 



(5? + 9ffi + v/121.9t - 54g2g2 + Slg^) ( -3N/fc , ^ , 
-^1.24 = , ^ . ' ^ '"3 (46) 

3l8v^g2 / Vl21^-54gf.^^l.j /3^.2 ^ /^^ 2 + 9 2 + ^m^^ _ 54^2 2 ^ gl^^A ^ 

V 3gj-9g3 + -y/121g5'-54gjg^+81g| V V ^ / / 

(g2 + 9g2 + V l21.9t-54g2g2 + 8lg4 ) (^^g^fc + y/ 3.^+3(llg?+9gi+7l21gt-54g;gi+81g|)^^ -^ 

>/l.25 = , , — T^2 

^121gf-54gfg|+81g| ' ' ' ' 



2 , o 2 / loi 4 r.1 2 2 I 01 4 A f /T , . / 9fc^+3(llgf +9gi- Vl21gt-54g;g|+8T^)y2 



"^2,24 — , r^3 

•^2,25 = 1 r^2 

and their symmetric counterparts. 

Here the J matrix show rotation between states 1-4 and the hght states, giving particle production and possible 
preheating. The reason that the SU(3) states do not rotate is that the 3 Q's have the same SU(2)-charge, and the 2 
diagonal SU(3) generators have removed the phases between them. 

In fact, changing the assignments such that the quarks have spht SU(2)-charges wiU change something, even the 
eigenvalues, but it wiU not change that J is nonzero and preheating is possible. 



VII. ONE FLAT DIRECTION - SUMMARY 

For the 2 flat directions mentioned as the most obvious candidate to be the inflaton in % QLD" and LLE", we 
find no preheating. The reason |12| found differently with a toy model direction of 3 superfields was that it was 
rather special to have 3 VEV-fields and only 1 broken generator. When only 1 generator was broken, only one phase 
difference was removed, and the second phase difference gave the preheating. However, for LLE'^ and U'^D'^D'^ 2 
diagonal generators are broken and there is no preheating due to the diagonal generators. We think, inspired by 
0, it makes sense to split the involved fields in those connected to VEV's by the diagonal generators (from here: 
Sector 1), and those connected to the VEV by the off-diagonal generators (from here: Sector 2). In this case, and 
we suspect in most others, the structure of Sector 2, is that the massive states are Higgses, and they all have the 
same eigenvalue. Therefore rotation does not have an effect (in fact, rotation does not make sense, since one cannot 
distinguish the eigenstates). In QLD'^ though, they have different eigenvalues - some fields connected to the VEV 
through SU{2), others through S'[/(3)c. However, for each field it is either or. Any difference from this, should be if 
a field is connected to 2 VEV's, one by a SU{2) and one by a SU{3)c generator. For Sector 1, it would take more 
than 3 fields (or less than 2 broken generators). This is what happens in QLQLQLE'^ . We can gauge away 4 phase 
differences, but this leaves 2 phase differences that can give the preheating. There could also be a mixing between 
sectors, if a field was connected to 1 VEV by a diagonal generator and to another by an off-diagonal one. However, 
this seems impossible for a single flat direction. 
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VIII. VD^D", LLE" SIMULTANIOUSLY 

The 2 directions first presented can co-exist. In fact, tliere is no reason wfiy they should not both get large VEV's 
0. It is not so easy to argue why there should be no preheating - since now we have 6 Sector 1 fields, and only 4 
diagonal generators to break. 



(47) 



where A is the relation between the absolute value of the VEV's involved. The Lagrangian reads 

3 



1=1 



(48) 



where for field (f>i D'^ = 
potential now looks like 



denotes the covariant derivative. The 



(49) 



To remove mixed kinetic terms we must reparametrise 



^c3 



V2 

(^ + ^3)e^(-+;fe) 



72 

(Av9 + ei2)e'^"^ + ^\ 

Cl5 + *66 



72 

Cl5 - 



n/2 

(A(^ + ei3)e'^"=+^^ 
(^(^ + $14)6'^"'+^^ 



(50) 



After verifying that the mixed derivatives have indeed been removed, we use coordinate derivatives, and find the 
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remaining kinetic terms (those not to second order) to be: (to zero'th order in ip) 



16 

E 

i=l 



10 



i=2 



1^2 • 2 
7?i0^i 



16 



4 / 2 



+ E 

i=12 

14 . 

;— 19 ^ 



2 

QO-2 



3 3 

+ 2*^^ 0=1^ + 2^^'^^^2^ + \/3v?f7i^'i + V3Aipd2(n- 
It seems the 2 directions do not "sec" each other. We find (after antisymmctrising) 
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Ex. 
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This looks as if the 2 parts are completely separated. The mass matrix for the physical excitations appears as 



(51) 



(52) 
























Vo 

with eigenvalues 
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(53) 



gl + gl + ^ y'(4 + 9^2)2^4 _ 6(_4 + 9^2)52^1 + g^i^ 

gl + gl - i^(4 + 9^2)2^4 _ g(_4 + <^A^)glgl + Qgl\ 



Ml 



Mil = Mh = iJ^gW 



Ml = Ml = M7 = Ml = Ml = 2gl(p^ 



(54) 



Mt. 



Mt: 



Mt 



0. 



Even tfiough this indeed looks like a mixing between the two directions, again J = (16 by 16 matrix) and there is no 
particle production and therefore no preheating. 



15 



IX. QLD" + LLE" - AN OVERLAPPING DIRECTION 



While LLE" + QLD"^ shows that more than 1 phase is not enough to secure preheating, it is clear that the 2 directions did 
not overlap. There are flat directions more intimately connected. One example of this is QLD'^ + LLE'^ - with one L-field 
in common. This is exiting, since here the flat directions cannot just have one phase each. Q and D"^ must be from different 
generations and have the same (or opposite, if you like) colour charge. It is easy to show that flatness is independent of phase, 
and that the common fleld shall have a VEV that is the square root of the sum of squares of the VEV's from the 2 directions. 



< i/e >= Vl + ylVe*^' (55) 
< ji >= (fie 
< T >= (fie ^ 

where A is the relation between the absolute values of the VEV's. The Lagrangian, covariant derivatives and the potential 
looks as before. 

To remove mixed kinetic terms, we must re-parameterise 



VT+A^ ^2{l + A^) J 

(f = {Aifi + ii)e^ ^ ^ Ml + A);\^ 2V6+8A2+6A4 

V2 
V2 

^c2 ^ ziii±!ilie<^+(-^-^)-^) (56) 
v2 



c3 _ -62+i|l3 



s = -pz e 

V2 



= (\/l + ^V + C3)e 



^2 

-^(&+iC9) , 64-iCl5 



yiTA^ ^2(1 + ^2) 



V^(l + A^) 



There are quite many kinetic terms now, but they include 
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and no cross terms. The U-matrix is (after antisymmetrising) 

(4 + (..l-)o +5(1 + .!")^ 



t^l,6 = — f^6,l = 1^2,6 = —[^6,2 
f/l,7 = —Ut^i = U2,7 = —[^7,2 



5^(1+^2) 

(2 + 3A2)((4 + 6A')j + 5(1 + A^)&) 



t/3,6 = — 1/6,3 

^3,7 = — f^7,3 



5VT0(1 + A2)V3 + 4^2 _^3A4 

4(l--l')i, +:)(l + -l-)a 
5v^(l + 

V2(-l + A^)(4(-l + A^)7 - 5(1 + A^)a) 



fA,7 = —Ur^4 = 1/5,7 = — [^7 



5^/5(l + ^12)73 + 4^2 _^3^4 

:^')7+5(l- 

i^(l + A2) 



rr _ n -rr - u - -(6 + 4^^)7 + 5(1 + 



6 



6(-l + A2)7 

(V8,9 = —U9,8 — — T7Z Tol h <T 



5VT0(1 + A2)V3 + 4^2 +3^4 

5(1+^2) 



^8,15 = — f/l6,8 = ?/9,14 = —Ul4, ^ 



(1+A2) 

0. 

This looks as if everything mixes - certainly it does not look as if there are 2 separable parts. The mass matrix is quite 
complicated, but the structure is like this 



\ Usx? J^sxs J 

where M is the mass matrix (with no 0-entries) for the previously mentioned sector 1, while D is a diagonal matrix which 
is for the previously mentioned sector 2. So, the sectors are clearly separated. The sector 1 part, has 3 very complicated 
eigenvalues, and zero is eigenvalue with multiplicity of 4. The sector 2 part has eigenvalues (entries) ordered after the ^-fields: 
(Si,Q,2A^gl,2A^gl,2A^gl,2A^gl,2{l + A^)gl, 2(1 + A^)gl)ip^ . Also, it is important to notice that M is time-dependant. All 
the elements involving 6 or 7 are time dependent. 



^7x7 07x 
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The elements (symmetry implied) are 
Mi,i = 4ff?,Mi,2 = -2g?,Mi,3 = -2\/rT^<??,Mi,4 = ^^,Mi,5 = ^^ 



2Agf 



M2.2 =91+ gl A/2,3 = Vl + A^g! - 52), M2,4 = A((-g?/3) + gl), ^ 

M3.3 = {l + A^){gl + gl),MiA = ~lAVTTA^{g! + 3gl),M3,5 = -^AVTTA^gl 
A\gi + 9g^2+l2gl) 2{gl - 6gl)A' AA\g! + 3gi) 

M4_4 = ,Al4,5 = ,Ai5^5 — 



Ml 



Ml,7 = 



9 ' ■ 9 ' ' 9 

2gf ((1 + A^) sin[g^^|±g^ - 2a] + ^TTA^ (-Asin[ ^'l,^+^f.>;' - 2a] + sin[i^g4^ + 2a])) 



(2(1 - A') sin[^tzi±^ _ 2a] + VT+A^ (a(3 + 2A') sin[i|±M^ - 2a] + (2 + 3A') sin['-^^^ + 2a]) ) 



2 „2 



M2,6 = - 



V3 + 7^2 + 7^4 + 3^6 

{gf - 52^) ((1 + sin[S^^i±^ - 2a] + VT+A^ (-Asin[l§±|4^ _ 2a] + sin['-^i0yL + 2a])) 



V5V1 + A^ 



_ jg'l - 9l) (2(A^ - 1) sin[^^^^ - 2a] - ^/T+W (a(3 + 2A^) sin[m^ _ 2a] + (2 + 3^^) sin[ig±|fll + 2a]) ) 
^^'"^ ~ ^/T0V3 + 742 _^ 7^4 + 3^6 



(P? + 52^) ((1 + A^) M^^m^ - + ^/^T:4^((-Asin[^(fii^ ^ 2a] + sin[f±M^ + 2a]))) 
■Afs.e = 



M3,7 = 

M4,6 = 

M4,7 = 

M5,6 = 

M57 = 



x/5 

(5? + ffi) (2(1 - A^) sin[^^^i±|^ _ 2a] + Vl + A^ (a(3 + 2^^) sin[i^y:M^ - 2a] + (2 + 3A^) sin[l(|±^ + 2a)) ) 



\/T0 V3 + 4^2 + 344 

A(g? + 3g|) ((1 + 4^) sin[g^^|±^ - 2a] + ^/^T:4^ (-Asinj^gij^ - 2a] + sin[^g±^ + 2a])) 

i^/Z^/TTA? 

A{gl + Zgl) (2(1 ^ A^) sin[^t(i)±^ - 2a] + ./TfW (yl(3 + 2A^) sin[i|±|4^ - 2a] + (2 + 3^^) sin[f±|4^ + 2a]) ) 



3VT0V3 + 7^2 + 7^4 + 3A6 
2^ff? ((1 + A^) sin[g^i^ - 2a] + ./T+A^ (-Asin[ig±M^ - 2a] + sin[f±^ + 2a])) 

3^^TTA2 



lAgl (-2(-l + A') sin[^^^i±^ - 2a] + VT+^ (a(3 + 2A^) sin[^§±|A^ _ 2a] + (2 + 34^) sin[ie±M^ + 2a])) 



3\/3 + 7^2 + + 3A6 
Af6,6, Ma^T, Mjj are much more complicated and are omitted here. 

X. QLD'^ + LLE" - PREHEATING 

In this case, where the mass matrix is time dependent, we must redo eg 1321 and we find 

C'^C = B'^A^{AB + AB) = -B'^UB + B^B (61) 

and therefore 

J = i f^c'^C-^ - ^ C' CV^ 

= Uv^i~B^UB)^-^{-B' UB)V^] +;-{V^B' B^- ^B' BV^] (62) 
= Ji + J2(with the obvious definition). 
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It has been shown numerically that Ji has the following structure - treating everything but ^i, a, 7 as constants - 
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(63) 



where NZ stands for a nonzero element. Since the columns of B were chosen such that the first 3 columns represent the massive 
states of sector 1, the next 4 the massless states in sector 1, and the last 8 columns represent the states of sector 2 - in the order 
in which the eigenvalues were mentioned, this represents particle production from the rotation between the 3 massive states 
and the massless states of sector 1, and particle production from the rotation between 2 of the massive and the massless states 
of sector 2. Reassuringly, there is no particle producing rotations between (indistinguishable) states of the same mass. Also it 
is easy to see that J2 will not alter this picture, since the last 8 columns of B axe constant. This make the last 8 columns of B 
(and therefore of B^B) zero and since multiplying by diagonal matrices cannot chance zero entries, it is clear that at least the 
last 8 columns of J are identical to the last 8 columns of Ji. Therefore the J-matrix is definitely nonzero, and there is particle 
production from the rotating eigenstates. It was shown in 12| that in general all phases will have nontrivial dynamics, which 
is necessary for the conclusion that J is nonzero. 



XI. THE METHOD OF FIELD COUNTING - AND ITS LIMITS 



A simpler approach for determining if preheating is possible is to count the fields, establishing the number of broken generators 
and thus the number of Goldstones and the number of Higgses and by subtraction finding the number of remaining, physical 
light degrees of freedom [^. This can be done sector by sector. 

LLE'^ breaks SU{2) x (7(1) completely. It has 6 fields (real) in sector 1, it breaks 2 diagonal generators and thus have 2 
Higgses and 2 Goldstones. This leaves 2 light degrees of freedom, which corresponds to the fiat direction. (It is clear that 
the sum of the phases cannot be gauged away.) However, one needs to argue why the 2 Higgses cannot rotate between each 
other - since they have different eigenvalues a possible rotation would be physical, and why the fiat direction stays out of this 
rotation. In this case the fiat direction - understood as 2 real fields, that is the direction itself, and the field combination that 
is orthogonal to it in all superfields individually - stays constant. So it is clearly not rotating. Sector 2 is very easy. There are 
4 fields, and 2 broken off-diagonal generators and therefore 2 Higgses and 2 Goldstones. Since the 2 Higgses have the same 
eigenvalue, there will surely not be particle production in this sector. 

U'^D'^D'^ breaks SU{?>)c x U(l) to U(1)new- It has 6 fields in sector 1, it breaks 2 diagonal generators and thus have 2 
Higgses and 2 Goldstones. This leaves 2 light degrees of freedom, which corresponds to the fiat direction. However, as before, 
one needs to argue why the 2 Higgses cannot rotate between each other. The fiat directions stay out of this rotation. Sector 2 
is again very easy. There are 12 fields, and 6 broken off-diagonal generators and therefore 6 Higgses and 6 Goldstones. Since 
the 2 Higgses have the same eigenvalue, there will surely not be particle production in this sector. 

QLQLQLE'^ breaks SU['i)c x SU{2) x 17(1) completely. It has 14 fields in sector 1, it breaks 4 diagonal generators and thus 
have 4 Higgses and 4 Goldstones. This leaves 6 light degrees of freedom, corresponding to the 1 fiat direction and 4 additional 
light degrees of freedom to which there can be rotations which give preheating. In Sector 2 there are 38 fields, and 8 broken 
off-diagonal generators and therefore 8 Higgses and 8 Goldstones. Since each field is exclusively connected to the VEV by 
SU{3) or SU{2) the Higgses cannot rotate between each other. 12 fields are completely decoupled (those of Q, differing in both 
SU{<i)c and Sf/(2)-charge from the VEV).Indeed there is rotation to some of the remaining 10 states, but it is hard to argue 
exactly why and to how many, without doing the full investigation. 

LLE" + VD^D" breaks 5'L''(3)c x SU(2) x U{1) completely. It has 12 fields in sector 1, it breaks 4 diagonal generators and 
thus have 4 Higgses and 4 Goldstones. This leaves 4 light degrees of freedom, corresponding to the 2 fiat directions. Again one 
needs to argue why the 4 Higgses cannot rotate between each other. The fiat direction clearly stays out of the rotation. In 
Sector 2 there are 16 fields, and 8 broken off-diagonal generators and therefore 8 Higgses and 8 Goldstones. Since each field is 
exclusively connected to the VEV by SU (3) or SU (2) the Higgses cannot rotate between each other. 

QLD'' + LLE" breaks S'(7(3)c x SU{2) x (7(1) to SU{2)a. It has 10 fields in sector 1 (one complex field in common), it 
breaks 3 diagonal generators and thus have 3 Higgses and 3 Goldstones. This leaves 4 light degrees of freedom, corresponding 
to the 2 fiat directions. Again one needs to argue why the 3 Higgses cannot rotate between each other and why in this case 
the light fields corresponding to the fiat directions does rotate with the Higgses. In Sector 2 there are 18 fields and 6 broken 
off-diagonal generators and therefore 6 Higgses and 6 Goldstones and 4 are completely decoupled. This leaves 2 light fields that 
can rotate. One can also argue, that the 4 down fields in Q and 4 strange fields in D'^ must represent the 4 color Higgses and 
4 color Goldstones with no particle production. 
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It seems clear, that while this counting is a nice tool to look for opportunities for preheating and to exclude preheating 
especially in sector 2, it is still necessary to do the full analysis in the unitary gauge to draw firm conclusions - at least when 
it comes to the role of the fields corresponding to the flat directions themselves. 



XII. SUMMARY AND CONCLUSION 



For the conclusions on one flat direction, see section IVIII For 2 flat directions, we have found that particle production is 
possible in QLD'^ + LLE'^ but not in U'^D'^D'^ + LLE'^. The difference seems to be the presence of a common field in the 
former case, but not in the latter. We have also found that it is necessary to transform to the unitary gauge after identifying 
the Goldstones, in order to make correct conclusions. Counting fields and broken generators can give hints to whether there is 
particle production or not, but it is not sufficient for firm conclusions. 

Finally, we shall stress that what we have shown is that there will be particle production in the QLD'^ + LLE"^ case. However, 
the statement that both directions are Ukely to get large VEV's [1] has not been investigated in this paper. Neither has the 
very recent claim that even if non-perturbative particle production happen, the main decay mode will still be perturbative pH ]. 
Also, whether the rotation of the flat directions are fast enough for this particle production to lead to preheating and thus not 
giving the effect of delayed thermalisation is outside the scope of the present paper. We presume that the situation is close to 
the situation in J^j - and that there will be very significant particle production. However, as stated in [T^, the effect of SUSY 
breaking terms in the Lagrangian has not been taken into account. 

But we can conclude that in order to determine the role of SUSY flat directions in (p)reheating, it is absolutely necessary to 
determine which flat directions get the large VEVs (only a limited number of the countless flat directions can get large VEVs 
at the same time) and if many directions get a large VEV a numerical study will probably be necessary to determine if the role 
of some additional flat directions can be ignored. 
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